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Abstract 

In the present article we have obtained new set of exact solutions of 
Einstein field equations for anisotropic fluid spheres by using the Her¬ 
rera et al. [I] algorithm. The anisotropic fluid solutions so obtained join 
continuously to Schwarzschild exterior solution across the pressure free 
boundary.lt is observed that most of the new anisotropic solutions are 
well behaved and utilized to construct the super-dense star models such 
as neutron star and pulsars. 

Keywords: anisotropic fluids; anisotropic factor; Einstein’s equations; Schwarzschild 
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1 Introduction 

The first ever exact solution of Einstein’s field equation for a compact object in 
static equilibrium was obtained by Schwarzschild in 1916. The static isotropic 
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and anisotropic exact solutions describing stellar-type configurations have con¬ 
tinuously attracted the interest of physicists Herrera et al. [HIIT]. Tolman [2] 
has proposed an easy way to solve Einstein’s field equations by introducing 
an additional equation necessary to give a determinate problem in the form 
of some ad hoc relation between the components of metric tensor. According 
to this methodology, Tolman [2] has obtained eight solutions of the field equa¬ 
tions, and his important approach still continues in obtaining the exact interior 
solutions of the gravitational field equations for fluid spheres. Buchdahl [3] 
proposed a famous bound on the mass radius ratio of relativistic fluid spheres 
is 2GM/c^r < 8/9, which is an important contribution in order to study the 
stability of the fluid spheres. Also Ivanov [3] has given the upper bound of 
the red shift for realistic anisotropic star models which cannot be exceed the 
values 3.842 provided the tangential pressure satisfies a strong energy condition 
(p > Pr + 2pt) and when the tangential pressure satisfies the dominated energy 
condition {p > pt)- Buchdahl [5] has also obtained a non-singular exact solution 
by choosing a particular choice of the mean density inside the star. 

The theoretical investigations of realistic fluid models indicate that stellar 
matter may be anisotropic at least in certain density ranges (p > 10^^ gm/cm^) 
(Ruderman and Canute m) and radial pressure may not be equal to the 
tangential pressure of stellar structure. The existence of a solid core due to 
presence of the anisotropy in the pressure was thought of by type-3A super¬ 
fluid (Kippenhahm and Weigert [7]), different form of phase transitions (Sokolov 
i) or by others physical phenomena. On the scale of galaxies, Binney and 
Tremaine have considered anisotropies in spherical galaxies, from a purely 
Newtonian point of view. The mixture of two gases (e.g. ionized hydrogen and 
electrons or monatomic hydrogen) can be described formally as an anisotropic 
fluid (Letelier [TU] and Bay in m) .The importance of equations of state for 
relativistic anisotropic fluid spheres have been investigated by generalizing the 
equation of hydrostatic equilibrium to include the effects of anisotropy (Bowers 
and Liang [H])- Their study shows that anisotropy may have non-negligible 
effects on parameters such as maximum equilibrium mass and surface red-shift. 

The relativistic anisotropic neutron star models at high densities by means of 
several simple assumptions showed that there is no limiting mass of neutron stars 
for arbitrary large anisotropy which is studied by Heintzmann and Hillebrandt 
m- However maximum mass of a neutron star still lies beyond 3-4 Mq. Also 
the solutions for an anisotropic fluid sphere with uniform density and variable 
density are studied by Maharaj and Maartens [14] and Gokhroo and Mehra [15] , 
respectively. Most the astronomical objects have variable density. Therefore, 
interior solutions of anisotropic fluid spheres with variable density are more 
realistic physically. 

Many workers have obtained different exact solutions for isotropic and anisotropic 
fluid spheres in different contexts (Delgaty and Lake [IB], Dev and Gleiser |T7| . 
Komathiraj and Maharaj [13], Thirukkanesh and Ragel m. Sunzu et al. |2D] . 
Harko and Mak [21], Mak and Harko [22], Ghaisi and Maharaj [23], Maurya 
and Gupta [24l|25j|2^, Feroze and Siddiqui [27], Pant et al. [28l|29l|3^, Bhara 
et al. m, Monowar et al. [32], Kalam et al. [33], Consenza et al. [34], Krori 
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[55] . Singh et al. [35], Patel and Mehta [37], Malaver [351133], Escupli et al. [30] . 
Herrera and Santos [4T1|42], Herrera et al. [43l l44] h 

The present paper consists nine sections, Section |T] contains introduction; 
Section [3] contains metric, its components and the field equations. Section [3] 
embodies the solutions of anisotropic fluid spheres in different contexts. Sec¬ 
tion [3] contains the expressions for density and pressure are mentioned for each 
fluid sphere.Section [5] consists the various physical conditions to be satisfied by 
the anisotropic fluid spheres. The analytical behavior of the solutions under 
the physical conditions (mentioned in Section [5]) are mentioned in the section |5| 
Section [7] describes the evaluation of arbitrary constants involved in the fluid 
solutions by means of the smooth joining of Schwarzschild metric at the pressure 
free interface r = a. The stability of models is proposed in the section [8] and 
finally section 0 includes the physical analysis of the solutions so obtained along 
with the concluding remarks. 


2 Metric, components and Field equations 

The line element of static spherical symmetric space time in the curvature co¬ 
ordinates a;* = (t,r, x,^) can be furnished as below, 

ds^ = B^{r)dt^ - 'tp-\r)dr^ - r^{dx^ + sin^ xd^^) (2.1) 

Einstein’s field equations given as 

\r51 ( 2 . 2 ) 

where, k = 

The components of the energy momentum tensor for spherically symmetric 
anisotropic fluid distribution is postulated in the form: 

P/ = (c^P + Pt)viV^ - ptSf + [pr - Pt)XiX ^, (2-3) 

where u* is four-velocity Bv'' = Sq, x* is the unit space like vector in the direction 
of radial vector, x* = VB6\, p is the energy density, Pr is the pressure in direc¬ 
tion of X* (normal pressure) and ptis the pressure orthogonal to Xi (transversal 
or tangential pressure). Suppose radial pressure is not equal to the tangential 
pressure i.e. pr ^ pt, otherwise if radial pressure is equal to transverse pressure 
i.e. Pr = Pt, it corresponds to isotropic or perfect fluid distribution. Let the mea¬ 
sure of anisotropy A = K{pt —pr) ,which is called the anisotropy factor (Herrera 
and Ponce de Leon [3^). The term 2{pt — Pr)I appears in the conservation 
equations = 0 (where, semi colon denotes the covariant derivative) which 
is representing a force due to anisotropic nature of the fluid. When pt > Pr, 
then the direction of force to be outward direction and inward when pt < Pr- 
However, if pt > Pr, then the force allows the construction of more compact 
object when using anisotropic fluid than when using isotropic fluid (Gokhroo 
and Mehra m)- 
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(2.4) 


In view of metric (EIJ, the Einstein field equations (12.21) give 

/ SttG \ 1 — tj) Ip' ( SttG \ 2B'ip ^p — 1 

y ^2 Pj - ^2 ^ ^7 ~ Br ^ r2 




1 

B r 


2ij 


B" 


B' 

rB 


= 2 A- 


where “dash” denote the derivative with respect to r. 


(2.5) 


3 Classes of solutions 


The field Eq. (12.51) has two dependent variables B{r) and ipij-), therefore Eq. (12.51) 
can admit infinity many solutions for different choices of B{r) and '>p(r) but all 
these solutions may or may not satisfy the physical conditions for the fluid 
spheres, for a given B{r), the Eq. (12.51) reduces to first order ordinary differ¬ 
ential equation in "(ppr). For its physically valid solution we will have to choose 
the metric potential B{r) such that B{Q) is non zero positive finite. This is a 
sufficient condition for a static fluid sphere to be regular at the centre. 

Let us take B{r) of the form 

B{r) = D{1 — cor^)''', where (n ^ 0), D is positive arbitrary constant. (3.1) 


Maurya and Gupta [521 US] have already obtained all possible anisotropic 
solutions of Einstein field equations for n is positive integer > 1 with cq < 0 
and —l<n<0, co>0 with different anisotropic factor A. 

But in present problem, we consider Bir) = D(1 — CQr'^)'' with n = —1, —2 
and —3 and anisotropy factor 


A = 


Aoco(cor^) " 

[1 - (n -f l)cor2]i+" 


; where, cq and Aq are positive constants. 


(3.2) 


Herrera et al [I] have proposed an algorithm for all possible spherically sym¬ 
metric anisotropic solutions of Einstein field equations. 

By using the Herrera et al [T] algorithm, the equations (|2.5p reduces in the 
form as: 


7 + 



y 


6 

r 




and 


y{r) 


B'{r) 1 

2B(r) r 


On integrating (|3.3a|) . we can obtain p) as: 


iP 


-2 


y{r){l + 

-3- dv H 


y2(^)g/[(4/r2i/(r))-|-2y(r)]iir 


(3.3a) 


(3.3b) 


(3.3c) 
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and 


(3.3d) 


y{r) = 


[1 - (n + l)cor^] 
r(l — Cor^) 


which further supply 




An 


= -20 


-d4> — 


= -20 




-d(j) 


+ for n = —1 (3.4a) 


and 


^/> = 


</>/" 


2/(n+l) 
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/ 


n —1 


d(() — 


/ 


n —1 


^ra+l^n(rt+3)/(n+l) ^ J ^2 g{n-l) / (n+l) 


-d(j) 


+ 0l. 


4>f 


2-n 


g2/(n+l) ’ 

for n ^ — 1. 
(3.4b) 


where, (j) = Cq?"^, / = (1 — ^) and g = [1 — (n + 1)^]. 

The Eq. (I3.4al) and Eq. (I3.4bl) give the following solutions: 

■tp = A.fPfe^'^ + (1 - 2(1)) f - (6 + Ao)f(l)e‘^^~‘^Ei{2 - 2(j>) + Aocpf; for n = -1, 

(3.5a) 


where Ei{2 - 2(j)) = log(2 - 2(j)) + wn 

N=1 

2 (15())^ - 25(()2 + 8) 15 


tP = 


)P = 


Ap){l + (p)"^ 
A. 




—(pf 

16 


2i(P^-(P + 2) + fil + P))^log 


16 
1 - 


(pil + (p) f log 


/ 


1 


Ap)f{l + 2(P)-^f 


P 


+ Ao(/>/^(l + 2^)(1 — 3(p + Scp"^) 
where h = (-81 + 130<() + 36(p^ - 720(P^ + 320<))4). 


1 + (^JJ 

320(;ilog 

; for n = —3; 


l + (p, 
for n = —2; (3.5b) 


(3.5c) 
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4 Expressions for energy density, pressures for 
different values of n 

(a) For n = —1 


+ 5(/) - 1) + (6 + Ao)/^(3 - 5(/) - - 20) 

+ (6 - 1102 203) + Ao[(-3 + 100 - 702) + 20e2-^-2/3K(2 - 20)]], 

(4.1a) 

= co[A/2(l + 0)e2^ + (2 - 100 + 702 - 203) - 6/2(1 + Ei{2 - 20) 

+Ao(l - 02)[1 - e^*-^fEi{2 - 20)]], (4.2a) 

(^Pt) = colAfil + 0)e2^ + (2 - 100 + 702 - 20^) - 6/2(1 + Ei{2 - 20) 

+Ao[(l + 0 - 02) - (1 - 02)e2^-2/i?*(2 - 20)]]. (4.3a) 

where, Ei{2 - 20) = - 2 E -. 

(b) For n = —2 




Co 


48 + 850 - 35002 9q^ 3 339^4 _ ^95^5 3O0/3(i + <^) 


- —f [45 - 150 - 28502-22503-Ao( 1503 + 1902 + 0- 3)] log ( -—^ 
16 ■' \1 + 0 

- ^[1103 - 1602 + 190 - 6 + 2/20(1 + 0)] 

(4.1b) 


Co 


+ Af (—3 + 0 + 190 + 150 ) 

1 


+(1 + 30)/3(1 + 0)2 

/3(l + 30)(l + 0)2log 


(15 + Ao) I I J.^2 


16 


(16 - 490 - 5002 90^3 30^4 _ 45^5) 

l-0\ 


l + ^J 


An 


+ T^/(1 + 30)(4 — 20 + 202) 
16 


(4.2b) 
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/^G 

V 


Pt = Co 


A{l + 3cj))f{l + ^f + 

, (15 + Ao),3 


(16 - 490 - 5002 9q^3 30^4 _ 45^5) 


16 


f{l + 30)(1 + 0)2 log (^Y^) 


+—^ (4 + 60 + 202 + 260^-60^) 
lo 


(4.3b) 


(c) For n = —3 


SttG 


= Co 


- A{3 - 150 + 9003 _ 405^4 447^5 _ 34 ^ 6 ) 


1 


— (21320 - 265650 - 5503402 455904^3 _ 44949^4 291200®) 

81 


- Ao/2(3 - 140 - 1002 + 9303 _ 99 ^ 4 ^ 


(4.1c) 


SttG 

-TrPr 


= Co 


Af{l + 70+ 1002)-(-356 + 4110+ 14602-382003-16800^ + 16000®) 


+ |^/‘(l + 2«(l + 5«log(10j 

+ Ao/2(l + 40 - 802 - 903 + 300^) 


(4.2c) 


SttG 


Pt = Co 


Af{l + 70+ 1002) + (356-4110-14602 + 382003 + 168O0'‘-16000®) 


QOf) 

+ + 50) log 


1 + 20 


243" ^ 

+ Ao {fil + 40 - 802 - 903 + 300^) + 03(1 + 20)2} 


(4.3c) 


5 Reality and Physical (well behaved) conditions 
for anisotropic solutions 

The physically meaningful anisotropic solution for the Einstein’s field equations 
must satisfy some physical conditions (Mak and Harko [5D], and Maurya and 
Gupta [22]): 

(i) The solution should be free from physical and geometrical singularities 
i.e. pressure and energy density should be finite at the centre and metric 
potential B{r) and tpir) have non zero positive values. 
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(ii) 


(iii) 


(iv) 


(v) 


(vi) 


(vii) 


(viii) 


(ix) 


(x) 


The radial pressure Pr must be vanishing but the tangential pressure pt 
may not vanish at the boundary r = ra of the fluid sphere. However, the 
radial pressure is equal to the tangential pressure at the centre of the fluid 
sphere. 

The density p and radial pressure pr and tangential pressure pt should be 
positive inside the star. 

[dpr/dr)r=o = 0 and [d’^Pr/dr^)r=o < 0 so that the pressure gradient 
dpr/dr is negative for 0 < r < r^. 

{dpt/dr)r=o = 0 and {d'^pt/dr'^)r=o < 0 so that the pressure gradient 
dpt/dr is negative for 0 < r < Va- 

{dp/dr)r=o = 0 and {d^p/dr'^)r=o < 0 so that the density gradient dp/dr 
is negative for 0 < r < r^. 

Conditions (iv)-(vi) imply that pressure and density should be maximum 
at the centre and monotonically decreasing towards the surface. 

Inside fluid ball the speed of sound should be less than that of light i.e. 


0 < 



< 1 , 0 < 



< 1 


In addition the velocity of sound monotonically decreasing away from the 
centre, the velocity of sound is increasing with the increase of density i.e. 


d 

dr 


{^) < > Oand f ($) < Oor (|j^) > 0 for 0 < r < r, 


In this contexts it is worth mentioning that the equation of state at ultra- 
high distribution has the property that the speed of sound is decreasing 
outwards (Canute [50]!. 


A physically reasonable energy-momentum tensor has to obey the energy 
conditions p> Pr + ‘^Pt (strong energy condition) and p + Pr + ‘^Pt > 0 . 


The red shift at center and at surface Za should be positive, finite and 
both bounded. 


The anisotropy factor A should be zero at the center and must be increas¬ 
ing towards the surface. 
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6 Physical properties of the new solutions 


(a) For n = —1 

The expression for pressures and density at the centre are as: 

Aoe 2 -K( 2)(6 +Ao) 


/SttG 


/SttG 


r—0 


(SttG 
= 1 — 


r—0 


— Cq [ A 2 


V c2 


r =0 


— Co [ —-A + 6 + 


i;i(2)(18 + 3Ao) - Se^Ao 


(6.1a) 

(6.2a) 


The pressure and density should be non zero positive at the centre and conse¬ 
quently A satisfy the following inequality: 


A>i-2 


Ei{2){6 + Aq) — Aoe^ 


and A < ( 6 -I- 


K(2)(18-f 3Ao) -3Aoe" 


(6.3a) 


The ratio of pressure-density should be positive and less than 1 at the centre 
i.e. < 1 which gives the following inequality, 

( pA - ( Pl\ - ( + 2) + Ape^ - £;z(2)(6-f Aq) \ ^ 

\pc2/r=0 Ve^(-^ + 6) + £^i(2)(18-f 3Ao)-3Aoe2y - 

(6.4a) 

By Eq. (I6.4al) . we get 

^<^^E.( 2)(6 + Ao)- 2 A„e^ 

Using the Eqs. (I6.3al) and (16.Sal) . A satisfies the following inequality: 


/ ^ ^ Ez(2)(6 + Ao) 


Ape" 


<A< 2 


Ez(2)(6-b Ap) 



Ap > 0. 
(6.6a) 


Differentiating (|4.2a|) with respect to r, we get an expression for the pressure 
gradient: 


/ 87rG dpr \ 
\ dr ) 


2 cpr 


Ae'^‘^{l -A(j) + (l)^ + 203) -f (-10 -f 140 - 603) 


- 6/3(1 -h 50)e3‘^-3Ef(2 - 20) 

- 6(5 - 120 - 303 -h lO03)e3'^-3£;f (2 - 20) 

-f Ap[-20 -f e^‘^-^Eii2 - 20)(-2 -f 50 - 30^) 


_g30-2y2(i_^0)^-(2_20)] 


(6.7a) 
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Thus it is found that extremum of Pr occurs at the centre i.e. 


pj, = 0 => r = 0 and 

^ ip'X=o = [(10 - + 6^*(2) + 30K(2) + Ao[2K(2) + Ei{2)]] 

^ ^ (6.8a) 


This shows that the expression of right hand side of equation (I6.8al) is negative 
for all values of A and Aq satisfying condition (I6.8ap . Then radial pressure (pr) 
is maximum at the centre and monotonically decreasing. 

Differentiating (I4.3al) with respect to r, we get 


/ 87rG dpt \ 
V dr ) 


= 2cor 


- 4(/) + + 2X + (-10 + Ucj) - + Ao(l - 2(j)) 


+ e^‘^-^Eii2 - 2(j))[Aoi-2 + 50 - 30^) - 6(5 
- (6 + Ao)f{l + 50)e2^-2K(2 - 20) 


120 - 30 ^ + 100 ^)] 
(6.9a) 


Which suggest that the extremum value of pt occurs at the centre i.e. 


Pf = 0 ^ r = 0 and (6.10a) 

^ {Pt)r=o = [(10 -A- Ao)e^ + &Ei{2) + 30K(2) + Ao[2K(2) + Ei{2)]] ; 

(6.11a) 


Under the condition (I6.6al) . the expression of right hand side of Eq. (IG.llal) is 
negative for all values of A and Ag. This shows that the tangential pressure (pt) 
is maximum at the centre and monotonically decreasing. 

Now differentiating equation (I4.1al) with respect to r , we get 


87rG dp 
dr 


2 cor 


Ae^^{5 - 230 ^ + 100 ^ + 80 ^) 


+ (-30 - 240 - 6003 _ 4804)e2d-2^j(2 - 20) 

+ 6/2(3 - 50 - 402)e2‘^-2K(2 - 20) - 220 + 60^ 
+ Ao[e‘^^-‘^f{Ei{2 - 20)(-5 - 190 + 120^ + 80^) 
- /(-5 + 90 + 802)Uf(2 - 20) 


- 2f(j)Ei{2 - 20)} - 2(-5 + 70)] ; 


(6.12a) 


where X{2 - 20) = ^ ^ 

The extremum of p occur at the centre i.e. 


p' = 0 => r = 0 and 

^ ip")r=o = [-(5^ + 10Ao)e" + Ei{2){30 + 5Ao) - K(2)(18 + 5Ao)] ; 

(6.13a) 
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Thus, the expression of right hand side of (I6.12al) is negative. This shows that 
the density p is maximum at the centre and monotonically decreasing towards 
the pressure free interface. 


The square of adiabatic sound speed at the centre, 




1 fdpr\ _ /(10-A)e2 + 6St(2) + 30K(2) + Ao[2K(2) + £;i(2)]\ 

? \ ^)r=o ~ V(5^ + 10Ao)e2 + K(2)(30 + 5Ao) + Ei{2)(l^ + 5Ao )) ’ 

(6.14a) 

1 fdpt\ Ao)e2 + 6K(2) + 30£;i(2) +Ao[2K(2) + K(2)]\ 

^ \ ^)r=o ~ V (5^ + 10Ao)e2 + Ei{2){m + 5Ao) + K(2)(-18 + 5Ao) ) ' 

(6.15a) 


The causality condition is obeyed at the centre for all values of constants under 
the condition (I6.6ap i.e. radial and tangential velocity of sound are monotoni¬ 
cally decreasing and less than 1. 


(b) For n = —2 

The central values of pressure and density are given by 


^ SirGpr ^ 


/ STrGpt \ 


Co 


r =0 




4 , — — (4A -1-8-1- Aq) 

/ r =0 ^ 




(-12A-h24-H3Ao) 


(6.1b) 

(6.2b) 


The central values of pressure and density should be non zero positive and 
finite, then A satisfies the following conditions:: 


A> -2- 


Ar 


and A <[ 2 — 


(6.3b) 


Subjecting the condition that positive value of ratio of pressure-density and 
should be less than 1 at the centre i.e. < 1 which leads to the following 

inequality. 


Pr 

pc^ 


r—0 


Pt 

pc^ 


r—Q 


By the inequality (I6.4b|) . we get 


1 + 


4A -|- 8 -)- Aq 
-12A-t24-t3Ao 


Ar 


< 1 


By the inequalities (I6.3bl) and (I6.5b|) . A satisfies the inequality 


Ar 


-2-^) <A< (1 + ^), Ao>0. 


(6.4b) 


(6.5b) 


(6.6b) 
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Differentiating Eq. (I4.2bl) with respect to r, we get 


SttG dpr 
dr 


= ^cir 


-2A. + log 


1 -, 


/ 2(-1 + 3<^+1302 + 9 (/) 3 ) 


64 + 130a; - 330x2 - 150x3 + 270a;4 ^ . 

---+ ^(2 - 160 + 1902 - 1003 - 304 ) 


Thus it is found that extremum value of Pr occur at the centre i.e. 


p'j. = 0 ^ r = 0 and 


SttG 


iPr)r=0 ~ 


-16^-64 + 2Ao 


(6.7b) 


(6.8b) 


Thus the expression of right hand side of the Eq. (j6.8bl) is negative for all values 
of constants A and Aq satisfying condition (I6.6bl) and it is showing that the 
pressure {pr) is maximum at the centre and monotonically decreasing. 
Differentiating Eq. (I4.3bl) with respect to r, we get 


87rG dpt 


= 2 cnr 


c4 dr ° 




1 - 


1 + . 


/2(-1 + 30+ 1302 + 903) 


64 + 130x - 330x2 _ 270 x‘‘ 


+ ^(2 + 4302 - 1003 - 30 ^) 


Thus it is found that the extremum value of pt occurs at the centre i.e. 


(6.9b) 


= 0 => r = 0 and 


87 rG 


(P ").=0 = 


l&A - 64 + 2Ao 


(6.10b) 


Thus the expression of right hand side of the Eq. (jb.lObI) is also negative for 
all values of A and Aq satisfying condition (j6.6bl) . This behavior shows that 
transversal pressure {pt) is maximum at the centre and monotonically decreas¬ 
ing. 

Now differentiating equation (I4.1b() with respect to r, the expression of den¬ 
sity gradient as: 


/ SttG dp 
\ c2 dr 


= 2 cnr 


55 


8(1 - 02 ) 


I r, A 15 + Aq 1 — < 
- { 2A-\ ---log ■ 


1 + . 


—(-30 + 880 - 7202 203 + 2503) 


/2(-5- 170 + 2502 + 4503) 
(6.11b) 


where, 5 = (32 + 860 + 1002 - 17603 - 3120^ + 9003 2730®). 

Thus the extremum value of p occurs at the centre i.e. 


p' = 0 =+ r = 0 and {p'')r=o = 2C( 


180 + 80A + 30Ao 


> 0; (6.12b) 
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The right hand side of Eq. (I6.12bl) is showing positive due to the inequality 
(|6.6bl) , and this condition gives that the energy density p is minimum at centre 
and monotonically increasing. 

and the square of its adiabatic sound speed at the centre, ^ (y^') ’ 

given by 


1 f dpr 


1 f dpt 


(16A-64 + 3Ao) 

c" V dp Jr=o ^ cH rfp y ~ [180 + 80A + SOAq] 


(6.13b) 


The causality condition is negative at the centre for all values of constants 
satisfying condition (I6.6bl) . Due to increasing nature of energy density, the 
solution is not well behaved for n = —2. 

Case 3: n = —3 

The central values of pressure and density are given by 


^ SirGpr ^ 


_ / STrGpt \ _ 


r—0 




J r—0 


Co [-^4 + 356 + A,o] 


/ 87rGp \ _ £o r_ 


\ c 2 / 

\ '^ / r=0 


81 


[-243^4 + 21320 - 243Ao] 


(6.1c) 

(6.2c) 


The central values of pressure and density should be non zero positive and 
Hnite. Then A satisfies the following conditions: 


/21 320 

A > (-356 - Ao) and A < f- Aq 


(6.3c) 


Subject to the condition that the ratio of pressure-density is positive and 

/hich leads to the i 

81 [A -h 356 Ao] 


less than 1 at the centre i.e. < 1 which leads to the following inequality, 




/ r=0 \P^ / r=0 

Eq. (|6.4c|) leads to 


-243A + 21320 - 243Ar 


< 1 


A < 


-f 


7516 


V 324 


By using the Eqs. (I6.3cl) and (I6.5c|) . we get the inequality for A as: 
- (356 -I- Ao) < A < - 

Differentiating (I4.2cl) with respect to r, we get 
87rG dpr 


(6.4c) 


(6.5c) 


(6.6c) 


dr 


= 2cnr 


Api{(j)) - (411 - 292(1) - lUmcf)^ - 6720(j)'^ + 8000(()^) 


320 


fO + 50 ) 


-^^P2(0) log + Ao(2 - 300 -I- 330^ -I- 16O0^-345a;'‘ -I-1800®) 

243 1 — (p 

(6.7c) 
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where, 


= (3 - 2A(t>-64>‘^ + 132(/.3-165(/)^ + 60(/)^), 
P2(</)) = (15 - AOcj) - 2520^ - 3360^ + 


Thus it is found that extremum value of Pr occurs at the centre i.e. 

p;=0^r = 0 and (p")r=o = 2cg (^3A - 411 + ^ + 2Ao^ . (6.8c) 


So the expression of right hand side of Eq. ()6.8c|) is negative for all values of 
A and Aq satisfying the condition ()6.6cp . This shows that the pressure {pr) is 
maximum at the centre and monotonically decreasing. 

Differentiating (I4.3cl) with respect to r, we get 


/ SnG dpt \ 
dr ) 


2cnr 


A(3 - 2A(l)-6(f>‘^ + 132(^3_j^g5^4 


-(411 - 292cl) - - 6720(/>3 + 8000((i'‘) + ^/^(l + ^c/)) 

+ ^(15 - 40^ - 252^2 _ 33g^3 j^20(/)^) log 

243 1 — d) 


+Ao(2 - 30<^ + 36((|2 + 176(()3-325a;‘‘ + 1800®) 


(6.9c) 


Thus it is found that extrema of pt occurs at the centre i.e. 
= 0 ^ r = 0 and^^ {Pt)r=o = ^ 


Thus it is clear that the expression of right hand side of Eq. (I6.10c|) is neg¬ 
ative for all values of A and Aq satisfying the condition (I6.6cl) . This condition 
gives that the transversal pressure is maximum at the centre and monotonically 
decreasing. 

Now differentiating equation (I4.1cl) with respect to r we get 


/87rG dp\ 

V d(j)) 


= 2cgr A(15 - 270^2 -f 660x® - 585a;'‘ -f 204a:®) -f ^ 

81 

f 320 

-|-(9 — 140a: -I- 540a:® — 495a:"^ -f 156x®) < log 

I 81 

-f Ao(2 - 300 -f 33a:® -f 16O0®-345x‘‘ -f 1800®) 


1 -f 20'( 

T^J 

(6.11c) 


pi = (-26885-55O180-f646860®-11497720®-H6856120^-29696O0®-f 1289600®) 
Thus the extrema of p occur at the centre i.e. 

p' = 0 ^ r = 0 and (7'")r=o = ^Cp ^ 26885 2Ao^ (6.12c) 
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Thus, the expression of right hand side of Eq. (I6.12cl) is negative for all values 
of A and Aq satisfying the condition (l 6 . 6 cl) .Then density p is maximum at the 
centre and monotonically decreasing. 


The square of adiabatic sound speed at the centre, 


^ (^)r=o’ gi^enby 


1 = 1 = / -99553 - 729A + 486Ao \ 

c^\dp)r=o V3[-26885 + 1215A+162Ao]y 

The causality condition is less than 1 and positive at the centre for all values 
of constants for all values of A and Aq satisfying condition (I 6 . 6 cl) . 


7 Boundary conditions for evaluation of constants 
A, D and cq 

The above system of equations is to be solved subject to the boundary condition 
that radial pressure Pr = 0 at r = a (where, r = a is the outer boundary of 
the fluid sphere). It is clear that m{r = o) = M is a constant and, in fact, the 
interior metric can be joined smoothly at the surface of spheres (r = a), 
to an exterior Schwarzschild metric whose mass is the same as in above i.e. 
m(r = a) = M (Masiner and Sharp [45]). 

Then the interior metric of this fluid spheres should be joined smoothly with 
Schwarzschild exterior metric such as B^{a) = 1 — 2u, where u = M/a, where 
M is the mass of the fluid sphere as measured by the exterior field and a is the 
boundary of the sphere. 

By joining (|3.5a|) and ()3.5c|) on the boundary of the anisotropic fluid spheres 
(r = a) and by setting (j>a = Coq'^ and ipanisia) = 1 - 2uams, we get the 
expressions of mass for n = — 1 and —3 as: 





- 3Aa + 24>l) 


+ (6 + Aoc/a)Ae^^^-^Eii2 - 2<^,) - Ao^af! 


(7.1) 


Ml 




+ + 2 <^a) log - M^Aoc/a.f! 

243 1 (pa 

Ml = (292 - SObc/a - 662^2 + 1796^^ - 1360()>^ + 320^®), 

M2 = il+2^a)il-^a+34'l). 

5 = (16 + - 5 oAa + ml + ml - ml)- 


(7.2) 


The arbitrary constant A in the expression (13.Sal) and (I3.5cl) can be determined 
by using the radial pressure Pr is zero at the boundary, the expressions of con- 
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stant ^ for n = — 1 and —3 are given as: 


A = 


[6(1 + 5M + Aofa] fy^^-^Ei(2 - 2</>a) - (2 - 10<(. + 7<tP - 2cj>^) - Aq/^ 


A = -^ 


/a(l + ^a)e 2 ‘^“ 


(-356 + 411^0 + imcjfi - m20(j)l - 1680-^4 1600(()5) 


(7.3) 


-|^/i(l + 2 ^.)(l + 5^.)logi^ 


- Ao/2(l + 40 - 802 - 903 + 3004) 


/4(1 + 70, + 1002) 


(7.4) 


Also the arbitrary constants D in the metric potential for the case n = — 1 
and —3 can be computed by the condition B{a) = 1 — 2uanis, 

D = fl[Af)afae^^^ + (1 - 20J - (6 + Ao)facl>ae^^^-^Ei{2 - 2cfa) + Ao0a]; 

(7.5) 


D = fa A0,/2(1 + 20,)-^/, 

+ Ao0a(l + 20a)(l — 30, 


^ - 3200a log 

J a 


30^) 


1 + ‘^4>a 


(1 + 20a) 


(7.6) 


The positive constant cq can be calculated by taking the surface density 2 x 1044 
gm/cm^ and using the condition cq = {^^a'^Pa) / (1 — 0a — a0a) for n = — 1 
and —3. 
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Figure 1: Variations of radial pressure Pr and the trace D — — 2Pt of the 

energy-momentum tensor for n = —1, c^af = 0.0829, Aq = 4.2395 and n = —3, 
coa2 = 0.0028, Aq = 0.7094 
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Tables for Numerical Values of physical quantities 

In Tables 1-2: Z = red shift, Solar mass Mq = 1.475 km, G = 6.673 x 
10“®cm^/gs^, c = 2.997 x 10^° cm/s, D = (8TrG/c^Co)p, Pr = (87rG/c‘^Co)pra^5 

P, = (8^G/cVo W, 7 = 



Figure 2: Variations of the radial velocity Vr = and ratio of radial 

pressure and density Rr = ^ of the energy-momentum tensor for n = — 1, 
coa^ = 0.0829, Aq = 4.2395 and n = -3, = 0.0028, Aq = 0.7094. 



Figure 3: Variations of the tangential velocity ft = y and ratio of tangen¬ 
tial pressure and density Rt = of the energy-momentum tensor for n = — 1, 
coa^ = 0.0829, Aq = 4.2395 and n = -3, = 0.0028, Aq = 0.7094. 
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r/a 


Figure 4: Variations of the anisotropy factor Ai=AxlO ^ and anisotropy 
factor A 2 = A X 10“® of the energy-momentum tensor for n = — 1, = 

0.0829, Ao = 4.2395 and n = -3, cqo^ = 0.0028, Aq = 0.7094. 


Table 1: n=-l, Ao= 4.2395, coa^=0.0829, Radius (a) = 16.0780 Km, Mass (M) 


= 1.7609 Me 


r/a 

Pr 

Pt 

D 

A 

/ dpj_ 
v c'^dp 

/ dpt 

V c^dp 

Pr 

Pt 


0.0 

1.3165 

1.3165 

20.7045 

0.0000 

0.9999 

0.9394 

0.0636 

0.0636 

0.2570 

0.2 

1.2543 

1.2683 

20.5763 

0.0140 

0.9918 

0.9312 

0.0610 

0.0616 

0.2554 

0.4 

1.0723 

1.1285 

20.1933 

0.0562 

0.9680 

0.9070 

0.0531 

0.0559 

0.2504 

0.6 

0.7852 

0.9117 

19.5603 

0.1265 

0.9292 

0.8672 

0.0401 

0.0466 

0.2421 

0.8 

0.4167 

0.6417 

18.6860 

0.2249 

0.8761 

0.8125 

0.0223 

0.0343 

0.2303 

1.0 

0.0000 

0.3515 

17.5839 

0.3515 

0.8092 

0.7427 

0.0000 

0.0200 

0.2151 



Figure 5: Variations of the red-shift {Z) of the energy-momentum tensor for 
n = -1, coa^ = 0.0829, Aq = 4.2395 and n = -3, coa^ = 0.0028, Aq = 0.7094. 
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Table 2: n=-3, Ao=0 .7094, coa^=0.0028, Radius (a) =3.1274Km, Mass (M) 
=0.8672Me _ 


r 

a 

Pr 

Pt 

D 

A X IQi® 

/ 

V c^dp 

V c-^dp 

Pr 

Wp 

Pt 

Z 

0.0 

4.0619 

4.0619 

13.1901 

0 

0.5160 

0.5160 

0.3079 

0.3079 

0.8597 

0.2 

3.8982 

3.8982 

13.1838 

0.00997 

0.5158 

0.5158 

0.2957 

0.2957 

0.8595 

0.4 

3.4077 

3.4077 

13.1650 

0.63900 

0.5153 

0.5153 

0.2588 

0.2588 

0.8590 

0.6 

2.5922 

2.5922 

13.1336 

7.2949 

0.5144 

0.5144 

0.1974 

0.1974 

0.8582 

0.8 

1.4549 

1.4549 

13.0896 

41.116 

0.5132 

0.5132 

0.1111 

0.1111 

0.8571 

1.0 

0.0000 

1.5748*10"® 

13.0331 

157.48 

0.5116 

0.5116 

0.0000 

1.208*10-9 

0.8556 


8 Stability of the stellar structure 

For physically acceptable model, one aspect that the velocity of sound should 
be within the range 0 < = {dp/c^dp) < 1 (Herrera [17] and Abreu et al. [H])- 

In present models, the expression for velocity of sound at the centre is given by 
equations (6.13a), (6.14a) and (6.13c). We plot the radial and transverse velocity 
of sound in Fig. 3 and conclude that all parameters satisfy the inequalities 
0 = '^ar = (dpr/c^dp) < 1 and 0 = v‘^^ = {dpt/c^dp) < 1 , everywhere inside 
the star models. From equations (6.13a), (6.14a) and (6.13c), we found that 
\v^t ~ '^arl < 1 at the centre and proved that velocity of sound in monotonically 
decreasing throughout inside the star. Also 0 < < 1 and 0 < < 1, 

therefore Ht-vlrl < 1. Now, to examine the stability of local anisotropic fluid 
distribution, Herrera’s [46] proposed the cracking (also known as overturning) 
concept which states that the region, in which radial speed of the sound is 
greater than transverse speed of the sound, is a potentially stable region. 

In our proposed models, the models are stable with the radius 16.0780 Km, 
Ao = 4.2395, coo^ = 0.0829 for n = —1 and radius 3.1274Km, Aq = 0.7094, 
cqo? = 0.0028 for n = —3. 

9 Physical analysis and conclusions 

In the present paper, the new set of anisotropic exact solution of Einstein’s field 
equations we have presented by taking the metric potential gn = (1 — cor^)” for 
n = -1,-2 and —3 and specific choice of anisotropic factor A which involves 
the anisotropic parameter Aq. The obtained solutions are utilized to contract 
the super dense star models with surface density 2 x 10^^ gm/cm^. It is observed 
that solutions are satisfying all reality and physical conditions (mention its in 
SectionjS]) for n = — 1 and —3. But the solution is not compatible for n = —2 due 
to increasing nature of its density (Section [HI Case (b)). The anisotropic fluid 
sphere possesses the maximum mass and corresponding radius are 1.7609^/0 
and 16.0780 Km for n = — 1, = 0.0829, Aq = 4.2395 and 0.8672Me and 

3.1274Km for n = —3, = 0.0028, Aq = 0.7094. The red shift for n = — 1 
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and —3 are monotonically decreasing towards the pressure free interface r = a 
and found that the red shift at the centre {Zq) and at surface (Za) are: (i) 
Zo = 0.2570 and = 0.2151 for n = -1, (ii) Zq = 0.8597 and = 0.8556 for 
n = —3 for both strong energy and dominated energy conditions. In our models, 
the red shift is also satisfied the upper bound limit for the realistic anisotropic 
star models (Ivanov [3]) and its behavior is represented by the Fig. [S] The 
Tables [1] and [2] shows the numerical values of physical parameters. Fig. 1 shows 
that the fluid spheres satisfies the strong energy condition. The behaviors of 
velocity and pressure density ratio are given by the Fig.[5]and[31 Fig.[3]represents 
the increasing nature of anisotropy factor for the fluids spheres. 
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